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ABSTRACT 



The measured proper motion of Fornax, expressed in the equatorial coordinate system, is 
(fj, a ,H$) = (47.6 ± 4.6, —36.0 ± 4.1) mas century -1 . This proper motion is a weighted mean of 
four independent measurements for three distinct fields. Each measurement uses a quasi-stellar 
object as a reference point. Removing the contribution of the motion of the Sun and of the Local 
Standard of Rest to the measured proper motion produces a Galactic rest-frame proper motion 
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of (/z° rf , /if rf ) = (24.4 ±4.6, — 14.3±4.1) mas century -1 . The implied space velocity with respect 
to the Galactic center has a radial component of V r = —31.8 ± 1.7 km s _1 and a tangential 
component of V t = 196 ± 29 km s _1 . Integrating the motion of Fornax in a realistic potential 
for the Milky Way produces orbital elements. The perigalacticon and apogalacticon are 118 (66, 
137) kpc and 152 (144, 242) kpc, respectively, where the values in the parentheses represent the 
95% confidence intervals derived from Monte Carlo experiments. The eccentricity of the orbit is 
0.13 (0.11, 0.38), and the orbital period is 3.2 (2.5, 4.6) Gyr. The orbit is retrograde and inclined 
by 101° (94°, 107°) to the Galactic plane. Fornax could be a member of a proposed "stream" 
of galaxies and globular clusters, however the membership of another proposed galaxy in the 
stream, Sculptor, has been previously ruled out. Fornax is in the Kroupa-Theis-Boily plane that 
contains eleven of the Galactic satellite galaxies, but its orbit will take it out of that plane. 

Subject headings: galaxies: dwarf spheroidal — galaxies: individual (Fornax) — astrometry: 
proper motion 



1. Introduction 

The Local Group is a dynamic environment. Galaxies swarm in the gravitational well of the cluster 
and move about each other. The Earth-bound observer had only a one-dimensional view of this plethora 
of motions until a few decades ago: measuring the radial components of velocities was "easy," measuring 
the tangential components was hard. Most of the difficulty lies in measuring the proper motion of a galaxy 
— which must be done with respect to cosmic standards of rest, such as high-redshift galaxies and quasi- 
stellar objects (QSOs). Since the size of a proper motion decreases with increasing distance, all else being 
equal, measuring motions of even the nearest galaxies required time baselines of several tens of years using 
images taken with the ground-based telescopes or just years when using those taken with the Hubble Space 
Telescope (HST). 

The analysis of the astrometric data is complex. The available ground-based images were acquired 
with a variety of telescopes and imagers and under different atmospheric conditions and air masses and, 
thus, they contain non-uniform distortions. Understanding and characterizing these distortions is essential 
because the expected measured proper motion is very small — on the order of a few tens of mas century -1 . 
Although space-based data do not contain distortions caused by the atmosphere, the analysis of such data 
is still complicated because of several factors. 1. Useful reference points are scarce. Even bright and 
compact background galaxies have larger positional uncertainties than star-like objects. The field of view 
of HST is small, thus, a typical image contains too few such galaxies to serve as useful standards of rest. 
Instead, a QSO with its star-like point-spread function (PSF) is the standard of rest of choice. Therefore, the 
availability of QSOs within the appropriate brightness range behind the target galaxy determines the number 
of independent images — fields. 2. Stars are scarce, ranging in number from a few tens to a few hundreds per 
field. This scarcity limits the ability to measure the PSF, to transform stellar coordinates between epochs, 
and to determine the average proper motion of the galaxy. 3. The images still contain distortions that are 
inherent to the detectors and optics of HST. The expected motion of stars with respect to the QSO is on 
the order of a few thousandth of a pixel in several years, therefore, a model for these distortions should be 
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accurate to a sub-pixel scale. None of the existing models has such accuracy. Despite these obstacles, clever 
observations and data analysis have allowed the list of nearby galaxies with a measured proper motion to 
grow steadily. 



In alphabetical order, the list now contains: Draco (Scholtz & Irwin 19931). Canis Maior (Dinescu et al. 


20051 


)). Carina (Piateketal. 2003). Fornax (Piatek et al. 20021 


Dinescu et al. 2004), LMC (Jones et al. 


1994; 


Kroupa et al. 199- 


1; KrouDa & Bastian 1997; Anguita et al. 


2000; Drake et al. 2001; Pedreros et al. 


2002; 


Kallivavalil et all 2006a). M33 (Brunthaler et al. 2005). Sagittarius (Dinescu et al. 2005a). Sculptor 


(Schweitzer et al. 1995; 


Piateketal. 2006), SMC (Irwin et al. 19961 Kallivavalil et al. 2006b). and Ursa 


Minor (Scholtz & Irwin 


19931 Schweitzer et al. 1997; Piatek et al. 2005). This sample is becoming large 



enough to reveal the three-dimensional nature of motions in the Local Group. Notable is the measurement 
for M33, the Triangulum Galaxy, the third most luminous ga l axy in the Local Group at a distance from 



the Milky Way comparable to that for M31. iBrunthaler et alJ (j2005l ) use radio observations with the Very 



Long Baseline Array of two groups of water masers located on opposite sides of the galaxy to derive a proper 
motion on the order of 50 /ias yr _1 . This is the smallest proper motion of a galaxy ever measured. The 
uncertainties in the model of the rotation of M33 and not measurement errors dominate the uncertainty of 
this proper motion. 

Just as the motions of planets around the Sun contain information about the mass distribution in and 
the formation of the Solar System, space velocities of galaxies contain such information for the Local Group. 
In addition, the space velocities give galactic orbits, and these can demonstrate if the galaxies evolved in 
seclusion or if they influenced each other through close encounters, collisions, or mergers. 

Dark matter continues to evade direct detection. Its gravitational interaction with the luminous matter 
and light is the key observational clue for its existence. On a galactic scale, the rotation curve of a spiral 
galaxy or the velocity dispersion profile of an elliptical galaxy are such key observations, whereas, on the 
scale of galaxy groups, it is the velocity dispersion of the member galaxies. Knowing all three components 
of the space velocities of the galaxies provides a check on the assumption of velocity isotropy in the mass 
estimates based on the velocity dispersion. It also offers additional ways of constraining the amount and 
distribution of dark matter. For example, for the Milky Way and its satellite dwarf galaxies one may: 1. 
Derive a lower limit for the mass of the Milky Way from the assumption that a satellite galaxy is bound 
or derive an estimate of the mass of the Milky Way out to the radius of the satellites from their velocity 
dispersion. 2. Derive a limit on the amount of dark matter in a satellite galaxy that approaches the Milky 
Way closely enough that the Galactic tidal force is significant compared to the self-gravity of the galaxy. 

The space velocities of galaxies can also test more general features of models for a cold dark matter 
cosmology by: 1. Comparing the distributions of the o bserved orbital elements with those predicted by nu- 
merical simulations of the formation of the Local Group (jGill et al H booilMoore et al.ll2006h. 2. Inve stigating 
possible close encounters and collisions among the satellite galaxies (|Zhaolll998t iDinescu et al.ll2004f ) . 



Two interesting ideas about the formation of the satellite galaxies of the Milky Wa y challenge the picture 

that t hese galaxies are surviving substructure, predicted by cold dark matter cosmology. Lynden-Bell fe Lynden-Belll 
([19951 ) propose that some Galactic satellite galaxies and globular clusters move along similar orbits — they 
are members of a "stream" formed when the tidal force of the Milky Way tore a progenitor satellite galaxy 
into fragments. The fragments continue to move together on the orbit of the prog enitor in the absence of colli- 
sions, significant orbital precession, or dynamical friction. This is a testable idea: Lvnden-Bell fe Lvnden-BeI3 
(j 19951 ) predict proper motions for the members of propose d streams. The othe r testable idea is based on 
the "planar" distribution of the Galactic dwarf galaxies. Kroupa et al. (2005) note that the 11 galaxies 
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nearest to the Milky Way lie within a disk whose thickness-to-radius ratio is < 0.15. iKroupa et all (|2005[) 
argue that such a distribution of satelli tes is very unlikely for t he s ubstructure predicted by cold dark matter 
cosmologies, though this is disputed by Libeskind et al. ( 2005 ) and lZentner et aD |2005h . If the orbits of the 
satellites are indeed within the disk, the geometry of the plane predicts the coordinates of their poles. Note 
that the galaxies do not have to share the same orbit — be in a stream — for their orbital poles to be the 
same, all that is required is that their orbits must be in the same plane. 

This article is one of a series that is attempting to better-define the kinematics of the Local Group and 
reports a measurement of the proper motion for the Fornax dSph galaxy. The meas urement deriv e s from 
a larger set of data than and, thus, replaces the preliminary measurement reported in iPiatek et al.l (|2002h . 
Section [2] presents Fornax, emphasizing those of its physical and structural properties that are relevant to 
our analysis. Section [3] describes the observations and data and Section |4] explains the process of deriving 
the proper motion making use of several key performance diagnostics. The next three sections, Sj5j fjH and 
f}71 present the measured proper motion and space velocity, calculate the implied Galactic orbit, and discuss 
the implications of the motion, respectively. The final section, ^|8l is a summary. 



2. Physical and Structural Properties of Fornax 



Among the Galactic dSphs, Fornax is one of the "giants." The luminosity of Fornax, given in Table 1, 
is comparable to that of Sagittarius, at least three times greater than that of the next most luminous dSph, 
Leo I, and about two orders of magnitude greater than that of the least luminous — Draco. Fornax and 
Sgr are the onl y dSphs with con firmed globular cluste rs. Reviews of research on dwarf galaxies in the 
Local Group by Mateo ( 1998h and van den Bergh ( 2000 ). and a study of structure of the Galactic dSphs by 
Irwin &: Hatzidimitrioul (1995) are convenient sources for a side-by-side comparison of the dSphs. 

Discovered by Sh apley (jShaplevHl938al fbh. Fornax is at the celestial location (a, S) = (02 h 39 m 53!l, -34°30'16.0") 
van den Berghl 1999 ) in the J2000.00 equatorial coordinate system. This location corresponds to {I, b) = 



(237.245°, —65.663°) in the Galactic coordinate system. The galaxy is only about 21° on the sky from the 
Sculptor dSph. 



Saviane et al.l (|2000h derived a distance modulus of (m — M)o = 20.70 ± 0.12 for Fornax using the 
magnitude of the tip of the red giant branch. The data consist of ground-based observations in the B, V, 
and / bands of four overlapping fields, each with an effective area of 10.7' x 10.7'. The implied heliocentric 
distance is 138 ± 8 kpc. Our study adopts this d istance when calculatin g distance-dependent quantities. An 
earlier measurement of the distance modulus by Buonanno et al. ( 19991 ) . using archival HST photometry in 
the F555W and F814W bands of cluster 4 in Fornax and its surroundings, produced a similar value. 

Several studies hav e investigated the structure of Fornax using wide-field imaging. A "benchmark" 
investigation is that by Irwin &: Hatzidimitrioul ( 1995h . who used star counts from a photographic plate 
taken with UK Schmidt telescope to derive the key structural parameters. Their field of view is about 
36 deg 2 . In a similar manner, the study also derives structural parameters for seven other dSphs. Table 1 
lists the values for the ellipticity, position angle, and core and tidal radii from llrwin &: Hatzidimitrioul (|1995l ). 
This contribution adopts these values where they are needed. 

An analysis of deeper images, albeit with smaller areal coverage, bv IWalcher et all |2003h produces a 
similar set of structural parameters as that in llrwin fc Hatzidimitrioul (|l995l ) . The images cover an area of 
8.5 deg 2 and reach a limiting magnitude of 23.5 in the V band. The resulting isopleth map shows a steeper 
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decline in the projected surface density on the south and east sides of the galaxy ve r sus th e n orth and west 
sides — a behavior noted previously by several studies, for example, iDemers et al.l (|l994T) or IStetson et al 



(119981) 



More recently, Coleman et al. ( 20051 ) surveyed Fornax producing images with a total areal coverage of 
10 deg 2 in the V and I bands and reaching a limiting magnitu de V = 20. The surface density of stars selected 
to lie along the giant branch shows that: 1) stars beyond the llrwin fc Hatzidimitrioul ( 19951) tidal radius are 



found to the north-west and south-east of the galaxy, which is the orientation of the minor a xis and 2) there 



i s an additional overdensity of stars in a shell about 1.3° to the north-west of the galaxy. IColeman et al 



(|2004l ) argued that the asymmetry in the surface density profile noted in the previous paragraph is caused 
by another shell about 1 7' to the south-east an d that this feature is dominated by a stellar population with 
an age of about 2 Gyrs. Coleman et al. ( 2005[ ) propose that these two shells and th e other extra-tidal stars 
are the signature of the infall of a gas-rich galaxy with My ~ —8. Deeper imaging by Olszewski et al.l (j2006h 



confirms that the inner shell contains an excess of stars formed in a burst about 1.4 Gyr ago. The location of 
these stars in the color- magnitude diagram indicates that they have [Fe/H]= —0.7, which is higher than that 
expected for a low-luminosity dwarf galaxy and implies that these stars formed from enriched gas originating 
from Fornax. The conflict between the evidence for infall — shells and extratidal stars — and the high 
metallicity of the inner shell may perhaps be resolved by star formation stimulated by a minor merger. 



Mateo et al] (|199ll ) measured radial velocities for 44 stars in two fields located in the direction of Fornax. 
One field coincides with the center of the galaxy and the other is about two core radii away from the center in 
the south-west direction. After excluding 12 likely non- member stars from the sample, the measured radial 
velocities imply a heliocentric systemic velocity of 53 ± 1.8 km s -1 . The sample shows no evidence for a 
significant rotation around the minor axis. The velocity dispersion is 9.9± 1.7 km s _1 in the central field and 
12.0 ± 2.8 km s _1 in the outer — thus, the velocity dispersion does not vary significantly with the projected 
radius. The study derives a mass-to-light ratio (M/Ly) of 12.3 ± 4.5 in solar units for the center of Fornax, 
but notes that M/Lys in the range betwe en 5.3 and 26 are al so acceptable given the uncertainties in the 
fitted structural parameters. More recently. I Walker et al.l (|2006j ) measured radial velocities for an additional 



176 st ars scattered throughout the galaxy and combined these measurements with those in iMateo et al 



([19911 ) in their analysis. The combined sample with the most stringent exclusion of non-members implies 
53.3±0.8 km s _1 for the heliocentric radial velocity, which our contribution adopts when calculating quantities 
that depend on this parameter. The global velocity dispersion is 11.6 ± 0.6 km s _1 ; the velocity dispersion 
calculated in radial bins does not vary with projected radius. 

Piatek et al.l (|2002 ) and Dinescu et al. I (|2004l) independently measured the proper motion of Fornax. 
The measurement bv lPiatek et al.1 (J2002) of (/i Q , /j,g) — (49 ± 13, —59 ± 13) mas century" 1 derives from HST 
images of three distinct fields using STIS or WFPC2. For two of these fields, the data are at two epochs 
separated by 1 - 2 years. F or the third field, th e data are at three epochs with a total time baseline of 2 



years. The measurement of iPiatek et al.l (|2002l ) is preliminary because it was based on partial data while 
the observing program was still ongoing. The measurement quoted in this contribution, whic h derives from 



the co mplete set of data, replaces the older preliminary measurement. The measurement by iDinescu et al 
([20041 ) of (fj, a ,fis) — (59 ± 16,-15 ± 16) mas century -1 derives from a combination of ground-based and 



HST imaging. The time bas eline between epoch s varies between 20 and 50 years. Some of the space-based 
data are the same as those in lPiatek et al.l (|2002l ). but they are such a tiny fraction of th e whole dataset that 



the measurement in IDinescu et al.l ([20041 ) can be considered independent from that in IPiatek et~al1 (|2002h 
or this article. 
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3. Observations and Data 

The data consist of images taken with HST in three distinct fields in the direction of Fornax. Each 
field contains at least one known QSO. Figure [1] shows the configuration of the fields — each represented by 
a square — on a 60' x 60' section of the sky. The center of Fornax is at (X, Y) — (0', 0'). North is up and 
East is to the left. The dashed ellipse in the figure represents the core. Its semi-major axis, ellipticity, and 
position angle are those in Table 1. One of the three fields — the smallest in size — is within the core and 
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Fig. 1. — A 60' x 60' section of the sky in the direction of Fornax. North is up and East is to the left. The 
dashed ellipse represents the measured core. The three squares depict the fields that this article studies. From 
North to South, they correspond to the fields FOR J0238 - 3443, FOR J0240 - 3434, and FOR J0240 - 3438. 
An arrow indicates the direction of the positive Y axis of the CCD. 

close to the minor axis. An arrow emanating from a field points in the direction of the positive Y axis in 
the CCD. The name of each field and the equatorial coordinates of its center are in the first three columns 
of Table 2. 
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Fig. 2. — A sample image from the epoch 1999 data for the FOR J0240 — 3434 field. It is an average of 
two images at one dither position with cosmic rays removed. Each image of the lensed QSO lies within a 
cross-hair. North is in the direction indicated by an arrow and east is 90° counter-clockwise from north. 



The most nor thern fie l d is F OR J0240 — 3434. It is approximately centered on i mage A of a lense d 
QSO discovered bv lTinnevi (|l995h and located a t (a, S) = (02 I M m 07;73, -34°34'19.8") (|Yentis et al.lll992h . 
Image B of the QSO is 6.1" away from image A. lTinnevl(|l995h measured the B-band magnitude and B — H 
color of image A to be 19.00 ± 0.05 and 1.42 ± 0.08, respectively, and of image B to be 19.77 ± 0.05 and 
1.92 ± 0.08, respectively. The redshift of the QSO measured from both images is 1.4. 

HST observed the FOR J0240- 3434 field on March 10, 1999; March 8, 2001; and March 8, 2003 using 
the Planetary Camera (PC2) of the Wide Field and Planetary Camera 2 (WFPC2) and the F606W filter 
(columns 4 — 6 of Table 2). The observations produced 18, 16, and 16 images, respectively, each with an 
exposure time of 160 s (column 7 of Table 2). The position angle of 156.6° for the Y axis of the CCD is the 
same to within 0.01° for all of the exposures. Figure [2] is a sample image from the 1999 epoch. 

The other two fields shown in Figure Q] are, in order of decreasing declination FOR J0240 — 3438 
and FOR J0 2 38 — 3443. Again, each field contains and is approximately centered on a known QSO. 
Tinnev et all |l997h provide information on these QSOs. The QSO in the FOR J0240 - 3438 field is at 
(a, 8) = (02 h 40 m 38!7, -34°38'58") (J2000.0), has a B-band magnitude of 20.2, and a redshift of 0.38. The 



QSO in the FOR J0238 - 3443 field is at (a, 5) = (02 
magnitude of 20.2, and a redshift of 2.0. 



1 38 m 43!8, -34°43'53") (J2000.0), has a B-band 



- 8 - 




500 
X [pixel] 



1000 



Fig. 3. — A sample image from the epoch 2000 data for the FOR J0240 — 3438 field. It is a sum of three 
images at one dither position with cosmic rays removed. The cross-hair indicates the QSO. 



HST observed the FOR J0240 - 3438 field on January 31, 2000; January 25, 2001; and January 29, 
2002. It observed the FOR J0238 - 3443 field on March 8, 2000; March 8, 2001; and March 9, 2003. For all 
of these observations, HST used the Space Telescope Imaging Spectrograph (STIS) with no filter (50CCD). 
There are 24 images per epoch of the FOR J0240 — 3438 field, each with an exposure time of 192 s, for a 
total of 72 images. The position angle of the Y axis for all of the images is 32.4° to within 0.05°. There 
are also 24 images per epoch in the FOR J0238 — 3443 field for a total of 72 images. The average exposure 
times at the three epochs are 193 s, 192 s, and 189 s, respectively, and the position angle of the Y axis is 
69.7° to within 0.06°. Table 2 summarizes the above information. Figures [3] and 0] are sample images from 
the earliest observations in the two fields. 



3.1. Charge Transfer Efficiency 

During the readout of a CCD, the electrical charge present in a pixel at the end of an exposure transfers 
from pixel to pixel towards the transfer register. For STIS, this motion is in the positive Y direction, or "up," 
whereas for PC2, it is in the negative, or "down," direction. Once in the transfer register, the charge moves 
along the X direction to the circuit that measures its amount. The transfer of the charge is imperfect: CCD 
pixels contain "charge traps" which both capture and release charge. Thus, the quantity of charge passing 
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Fig. 4. — A sample image from the epoch 2000 data for the FOR J0238 — 3443 field. It is a sum of three 
images at one dither position with cosmic rays removed. The cross-hair indicates the QSO. 



through a pixel can either increase, if more charge is captured than released, or decrease, if the reverse is true. 
Because the charge traps have a finite capacity, they have a smaller fractional effect on bright sources than 
on faint sources. The origin of charge traps is not completely understood, but it is known that the efficiency 
of charge transfer in STIS and WFPC2 decreases with time ( Whitmore et al. 1999t Goudfrooii et al. 2006) 



most likely because of defects in the crystal lattice of the CCD created by cosmic radiation (jJanesick et al 



19911 ). Therefore, the effect of charge traps on the photometry and astrometry increases with time. 



Bristow et al.1 (|2005f ) describe the various effects caused by charge traps and why these effects arise, 
emphasizing the impact on time-dependent astrometry such as the measurement of a proper motion. Figured] 
illustrates how charge traps can create a spurious proper motion. Consider a field that contains three stars 
with similar signal values and a QSO with a much larger signal. In the figure, the stars and the QSO are 
represented by open circles and a solid square, respectively, and are labeled as 51, 52, 53, and Q. The 
QSO is the reference point. Assume that the stars, the QSO, and the observer are at rest with respect to 
each other, which means that the observer should detect no motion for any star with respect to the QSO. 
Figure [5] depicts two hypothetical images of this field taken at t = and at t = T with exactly the same 
pointings of the telescope and with a CCD that contains charge traps whose number has increased with 
time. Let the measured centroid of a source be depicted by the location of the source in an image. Because 
in the time interval T separating the two images the number of charge traps has increased, a centroid in 
the later image has a smaller value of Y — it has shifted away from the transfer register. Here we ignore 
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Fig. 5. — Hypothetical images taken at t = (left) and ait — T (right) with a CCD containing charge traps 
whose number increases with time. There are four sources in an image: three stars and a QSO. The open 
circles, labeled 51, S2, and 53, represent the three stars while the solid square, labeled Q, represents the 
QSO. The location depicted for a source in an image corresponds to its measured centroid. The signal in 
the QSO is much larger than the signal in each star, the latter assumed to be about the same. Because the 
number of charge traps has increased between t = and t = T, a measured centroid in the "later" image 
has a smaller value of Y, that is, the centroid has shifted in the direction of decreasing Y values. Arrows 
in the figure represent these shifts. The magnitude of the shift depends on the location of a source in the 
image and on its signal. 
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a possible shift in the X directio n resulting from ine fficient charge transfer in the serial register since this 
inefficien cy is negligible for STIS (jBrown et al.l 120021 ) and one-third as large as in the parallel direction for 
WFPC2 (jHever et al.l 120041 ). The arrows SyQ,Sysi,5ys2, and Sys3 represent these shifts for the QSO and 
the three stars, respectively. The magnitude of the shift is larger the smaller the signal of the object and the 
greater its distance from the serial register. Thus, for example, S3 has a spurious motion in the coordinate 
system of the CCD of (Sys3 — (5j/q)/T, which implies that the observer would measure a proper motion for 
S3 with respect to the QSO. 

Both WFPC2 and STIS contain charge traps whose numb er increases with tim e as a consequence of the 
ongoing production of lattice defects by the cosmic radiation (jJanesick et al.lll99ll ). The transfer register is 
at the bottom of the image for PC2 and at the top for STIS. Therefore, the centroids shift towards larger 
values of Y, or "up," for PC2 and towards smaller values of Y, or "down," for STIS. Accurate astrometry 
requires either that the shifts are accounted for in the analysis or th at a CCD image b e "rest ored" so that 
its pixels have the same values as immediately before the readout. Bristow fc Alexov ( 20021 ) developed a 
software package that restores an image taken with STIS, which was used to correct all of our images taken 
with STIS prior to our analysis. No corresponding software exists for restoring an image taken with PC2. 
In this case, our analysis measures and removes the shifts — the details are in Section I5TT1 



To estimate how greatly the decreasing charge transfer efficiency can affect the measured proper motion, 
this study performs a parallel analysis ignoring the effects of charge traps and comments on the difference. 



4. Measuring Proper Motion 



Piatek et al.l (|2006l ) and references therein describe our method of deriving a proper motion from un- 
dersampled images taken with HST. This article outlines briefly the four major components of the method. 

1. Determine the initial estimates for the centroids of objects — stars and the QSOs. This task uses 
stand-alone astronomy software packages such as DAOPHOT, ALLSTAR, and DAOMASTER |stetsonlll987 . 
1992LI1994I ). For each held and epoch, this step produces a list of objects at that epoch and the translational 
offsets between the chronologically hrst image of an epoch — the fiducial image — and those that follow. 
The centroid of an object on the list is in the coordinate system of the hducial image; however, its centroid 
in any other image derives from the known offsets. 



2. Construct an effective point-spread function (ePSF; Anderson fc King (j200Clh ) for each field and 
epoch. Using a select set of stars and the QSOs, our method builds a single, position-independent or constant 
ePSF that is used for all of the images of a field at an epoch. An object contributes to the construction of 
the ePSF if its signal-to- noise ratio (S/N) is greater than some limit — 15 for these data — and if its inner 
3x3 array of pixels (see iPiatek et al.l (|2002) for details about the structure of the data) have data quality 
flags equal to 0. If there are N images for a given field and epoch, an object may contribute up to N times 
to the construction of the ePSF — fewer if the object is flagged out in some images. 

3. Derive accurate centroids for objects by fitting an ePSF to the image of each object using least 
squares. The fitting procedure for a set of images at an epoch involves an interplay between the centroids, 
the ePSF, and the transformation between the coordinate system of an image and the fiducial image. An 
improvement in the centroids leads to more accurate transformations and a better estimate of the ePSF 
and this is iterated until a stable solution ensues. Averaging the centroids from each image transformed 
to the coordinate system of the fiducial image produces the coordinate for an object at one epoch and its 
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uncertainty. The fiducial coordinate systems for different epochs are generally offset from each other because 
of small changes in the telescope pointing and rotation and are also not congruent because of changes in 
image scale and telescope motion (the aberration of starlight). 

4. Determine the motion of the QSO in a "standard" coordinate system. In this coordinate system, 
the stars of the dSph are at rest while all other objects move. A proper motion of the dSph derives from 
the reflex motion of the QSO in this coordinate system. The standard coordinate system is defined as 
the fiducial coordinate system of the first epoch. The X and Y coordinates of a star of the dSph in this 
system are the average of its coordinates at the three epochs transformed into the standard coordinate 
system. A six-parameter linear transformation converts coordinates to the standard coordinate system. The 
transformation also includes a fitted linear motion for the QSO and, iteratively, for objects whose large 
scatter in their transformed coordinates indicate that they are not members of the dSph. For data taken 
with PC2, an additional term in the transformation that is linear with y and depends on the flux of an object 
corrects the effects caused by charge traps. 

Both STIS and WFPC2 exhibit geometric distortions which are several pixels at the edge of the field. 
These are large compared to the approximately 0.001 pixel precision in measuring the positions of objects 
required to determine the proper motion of a dSph. Even the accuracy with which the distortions are known 
is larger than the required precision, so our observations must minimize the impact of the distortions. Thus, 
the observations attempted to place objects on the same pixels of the CCD at each epoch. If the execution 
of this plan were perfect, no correction for the distortions would have been necessary. Although the roll 
angles are similar to a high degree, there are linear offsets between epochs on the order of a few tens of 
pixels. Therefore, o ur method corrects the distortions usin g the prescriptions f rom t he STIS Data Handbook 
(jBrown et al.l l2002) for the data taken with STIS and from I Anderson fc Kind (|2003T ) for the data taken with 
WFPC2. 



An error in the manufacture of t he WFPC2 CCDs makes every 34th row in the PC2 images narrower than 
the other rows (jShaklan et al.lll995l ). Dubbed the "34th-row defect," it can produce a spurious contribution 
to the measurement of the proper motion. Our me thod corrects the centroids for the effect of this defect 
using the prescription from Anderson fc King ( 19991 ). 



4.1. Flux Residuals 



Piatek et all (|2002h developed diagnostics for the performance of their algorithm that determines the 



centroids of objects. The remainder of this section describes the "flux residual" and "position residual" 
diagnostics and how the algorithm performed in the three fields. 

The flux residual diagnostic, 1ZJ 7 , defined by Equation 22 in IPiatek et alJ ( 2002h is a measure of how 



well the constructed ePSF matches an image of an object. If the match is perfect, 1ZJ- — 0. However, if the 
image of an object is wider than the ePSF, IHF > 0; otherwise, TZJ- < 0. If TZJ- does depart from zero, the 
departure increases with increasing S/N for an object. 

Plotting 1Z T as a function of position in a CCD can reveal potential flaws in the data and analysis. 
In the presence of only random noise, a plot of 1ZT versus X or Y would show that the points scatter 
around 1ZJ- = 0. Any other distribution is undesirable and can mean, for example, that the true PSF varies 
with position in a CCD, that the constructed ePSF has an incorrect shape, or that the PSF of an object is 
distinctly different from the ePSF, as it would be in the case of a galaxy with a resolved core. 
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4.1.1. Flux Residual Diagnostic for the FOR J0240 - 3434 Field 

Figure [5] plots 1ZJ 7 versus X in the left column and 1ZJ- versus Y in the right column for the FOR J0240— 
3434 field imaged with PC2. From top to bottom, the plots are for the epochs 1999, 2001, and 2003. Solid 
squares mark the points corresponding to image A of the QSO, whereas the solid triangles do the same for 
image B. All of the plots show a modest "bow-shaped" dependence of 1Z T on X and Y indicating that the 
true PSF varies with location. The constructed ePSF is narrower than the true PSF around the edges of the 
CCD and wider near the center. The values of H.T for both images of the QSO tend to be negative; those 
for image A are more negative than those for image B, most likely because of the variation of the PSF across 
the image. The 1ZJ 7 for the QSO is more negative than those of most nearby stars, which could be due to 
the higher S/N of the QSO or, perhaps, to a narrower PSF caused by the bluer color of the QSO. 

4.1.2. Flux Residual Diagnostic for the FOR J0240 - 3438 Field 

Figure [7J plots 1Z T versus X (left column) and 1ZJ- versus Y (right column) for the FOR J0240 — 3438 
field imaged with STIS. From top to bottom, the rows of plots are for the epochs 2000, 2001, and 2002. 
The values of 1ZJ 7 for the QSO are positive for all epochs and most of these points are outside of the top 
boundary in the upper two rows of plots. The PSF of the QSO is thus wider than the ePSF, most likely 
because, with a redshift of 0.38, this QSO is near enough for its core to be resolved. 

The values of TLT for stars in the top and bottom rows of plots are well behaved. They show no clear 
trends with either X or Y and scatter around 1ZJ- — 0. A few bright stars have TZTs comparable to those 
for the QSO. However, the plots in the middle row are not well behaved. The plot of 1ZJ- versus X shows a 
clear trend: IZFs tend to be positive for X < 500 pixel and negative otherwise. A line with negative slope 
crossing TZT = around X w 500 pixel could fit the trend. Although the plot of 7ZJ- versus Y does not 
show trends, the scatter around 1ZJ- = is larger than that in the corresponding plots at the other two 
epochs. The likely cause for this behavior is a dependence of the PSF with X. Remembering that the ePSF 
is a global average of individual PSFs, the PSF is wider than average for X < 500 pixel and narrower than 
average for greater values of X. A variation of the PSF with X in STIS seems isolated to this field and 
epoch. We are unable to trace the source of this dependence. Position residual diagnostics, 1ZX and IZy, 
discussed in Section [4.21 will show if the variable PSF had a significant impact on the measured centroids. 

4.1.3. Flux Residual Diagnostic for the FOR J0238 - 3443 Field 

Figure [5] for the FOR J0238 — 3443 field is analogous to Figure [7J From top to bottom, the rows of plots 
are for the epochs 2000, 2001, and 2003. 

None of the plots shows a trend between TZT and either X or Y . The values of TZ!F for the QSO tend to 
be negative — most distinctly in the bottom row of plots — indicating that the PSF of the QSO is narrower 
than the ePSF. 
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Fig. 6. — Left column: Flux residual 1ZT as a function of X. Right column: Flux residual 1ZT as a function 
of y. The plots arc for the FOR J0240 - 3434 field that was imaged with PC2. From top to bottom, the 
rows of plots are for 18 images at epoch 1999, 16 at 2001, and 16 at 2003, respectively. The figure shows 
points only for those objects that have S/N > 15. The solid squares mark points corresponding to image A 
of the QSO, while the solid triangles do the same for image B. 
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Fig. 7. — Left column: Flux residual 1ZJ 7 as a function of X. Right column: Flux residual 1ZJ 7 as a function 
of Y. The plots arc for the FOR J0240 - 3438 field, which was imaged with STIS. From top to bottom, the 
rows of plots are for 24 images at each of epoch 2000, 2001, and 2002, respectively. The figure shows points 
only for those objects that have S/N > 15. The solid squares mark points corresponding to the QSO. 
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Fig. 8. — Left column: Flux residual 1ZT as a function of X. Right column: Flux residual 1ZT as a function 
of Y. The plots arc for the FOR J0238 - 3443 field, which was imaged with STIS. From top to bottom, the 
rows of plots are for 24 images at each of the epochs 2000, 2001, and 2003, respectively. The figure shows 
points only for those objects that have S/N > 15. The solid squares mark points corresponding to the QSO. 
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4.2. Position Residuals 

The position residual diagnostic consists of plots of the position residuals 1ZX = (Xq) — Xq and IZy = 
(Yq) — Yq versus the pixel phase & x = Xq — Int(Xo) and Q y = Yq — Int(Yo)- Here (Xq) and (Yq) are the 
components of the mean centroid in the fiducial coordinate system — the system of the first image in time at 
a given epoch. The function Int(x) returns the integer part of a variable x. The coordinate transformation 
to a fiducial system involves translation, rotation, and scale change. The coefficients of the transformation 
result from a least-squares fit using objects common to the fiducial image and the one transformed to it. 

In the presence of only random errors, position residuals would scatter around 1Z X = or IZy = with 
constant rms for all values of $ x or <& y . Trends and structures in the plots are symptomatic of systematic 
errors. 



4-2.1. Position Residual Diagnostic for the FOR J0240 - 3434 Field 

Figure [9] shows plots of the position residuals versus pixel phase for the FOR J0240— 3434 field. Panels 
a, b, and c are for epochs 1999, 2001, and 2002, respectively. Only two plots, 1ZX versus $ x and TZy versus 

in Figure [9ji show evidence for trends among the points for images A and B of the QSO. In the first of 
these two plots, the values of TZX are mostly positive up to Q x ~ 0.8 and then negative for larger values of 
$ x . In the second, the IZys show a sinusoidal dependence on Q y . 

No plot for any epoch shows trends or structures for points corresponding to stars, though wc caution 
that the small number of points and their large scatter due to the modest S/N of the majority of stars would 
make such trends and structures hard to see. 

Even though Figure [6] implies that the PSF varies with position in the image at all three epochs for 
this field, Figure [9] shows that the QSO has pixel phase errors only at the first epoch. Thus, the origin of 
the pixel phase errors for the QSO are not clearly caused by the variation of the PSF with position in the 
field. Neither are they clearly caused by the variation of the PSF between images at one epoch because, as 
Figure [5] shows, that variation is similar at all three epochs. 

4.2.2. Position Residual Diagnostic for the FOR J0240 - 3438 Field 

Figure HO] plots TLX and TZy versus $ x and $ y for the FOR J0240 - 3438 field. Panels a, 6, and c are 
for epochs 2000, 2001, and 2002, respectively. 

None of the plots in Figure [TUb shows convincing evidence for systematic errors. However, the plots of 
TZX versus $ x in Figures [TUb and [TUb do. In the first plot, TZX has a linear trend with & x with positive 
slope. In the second plot, the trend is more sinusoidal. In both plots, IZXs corresponding to stars show 
greater scatter than IZys for the same objects. To a lesser degree, the plot of IZy versus $ y in Figure fTUh 
shows a sinusoidal trend too. 

As is the case for the previous field, there is not a clear correlation between a variation of the PSF with 
position and the presence of pixel phase errors for the QSO. For the 2001 epoch, both a variation of the 
PSF with X and pixel phase errors in X are present. In contrast, for the 2002 epoch, pixel phase errors are 
present even though the PSF does not vary with position. The 2002 epoch has better agreement than the 
2000 and 2001 epochs between the PSF of the QSO and that of the stars and it also shows less variation of 
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Fig. 9.— Plots for the FOR J0240 - 3434 field of the position residuals, TZX and Tiy, as a function of the 
pixel phase, $> x and <fr y . Panels a, 6, and c are for the epochs 1999, 2001, and 2002, respectively. Solid 
squares and triangles correspond to images A and B of the QSO, respectively. The plots show points only 
for those objects with S/N > 15. 



the PSF between the images of an epoch. 



4.2.3. Position Residual Diagnostic for the FOR J0238 - 3443 Field 

Figure ITU plots TLX and Tiy versus $ x and <$> v for the FOR J0240 - 3438 field. Panels a, b, and c are 
for epochs 2000, 2001, and 2003, respectively. 

No plot in the three panels shows clear evidence for systematic errors. The absence of systematic trends 
in these plots is consistent with the absence of trends in Figure [8] 



5. The Proper Motion of Fornax 



At this point in the analysis, t here are three lists of coordinates — one per epoch — for each field. 
Equations (1) - (4) in IPiatek et al. |2006h give the form of the transformation between epochs. This most- 
general linear transformation has six free parameters ci — Ce- The t ransformation also includes an additional 



term that is linear with Y (see Equations (7) and (8) in lPiatek et al.l (|2005l )) if the fitting procedure corrects 
for the effects caused by decreasing charge transfer efficiency. The slope of the linear correction, b, is an 
adjustable parameter and this correction is applied only to objects whose S/N is below a specified limit. 

The systematic errors in the object coordinates reflected in the trends in 1ZJ 7 , 1ZX, and IZy with X, Y, 
<f> x , and <& y will contribute to the error in the centroid of an object transformed to the standard coordinate 
system. Because we are unable to eliminate the source of these trends, our procedure propagates these 
errors into the uncertainty in the proper motion by increasing the uncertainties in the coordinates of all 
objects until the y 2 of the s catter aro und the transformat ion is equal to one per degree of freedom (see the 
discussions in lPiatek et al. ( 2005 ) and Piatek et al. ( 2006 )). 



5.1. Motion of the QSOs in the FOR J0240 3434 field 

The PC was the imaging detector for this field; therefore, correcting for the effects caused by decreasing 
charge transfer efficiency requir es fitting for the par ameter b. The coordinates of objects with a S/N < 30 are 



corrected using Equation (8) of IPiatek et al.l (2005). For objects with a larger S/N, the size of the correction 
decreases linearly with increasing S/N until a value of 100; beyond this value there is no correction. This 
field contains a lensed QSO with two images, A and B, and so it provides two measurements of the proper 
motion. 

Of the 233, 209, and 231 objects measured at the three epochs, respectively, 191 are common to all 
epochs. The greatest difference in the pointing of HST is between the 1999 and 2001 epochs, where the 
difference is about 20.6 pixel in the X direction and about 5.6 pixel in the Y . The choice for the individual 
X 2 that triggers fitting for uniform linear motion is 15. The multiplicative constant that produces a x 2 of 
one per degree of freedom is 1.304. The fitted value of b is —6.5 x 10 -5 . 



Figure [T2l plots the position residuals, RX and RY, defined by Equations (10) and (11) in lPiatek et al 



(|2005l ) as a function of position in the standard coordinate system. A residual is the difference between the 
average centroid in the standard coordinate system and the centroid measured at an epoch, corrected for any 
fitted motion, transformed into the standard coordinate system. The first digit of the subscripts on RX and 
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Fig. 10.— Plots for the FOR J0240 - 3438 field of the position residuals, TLX and Tiy, as a function of the 
pixel phase, & x and Q y . The panels a, 6, and c are for the epochs 2000, 2001, and 2002, respectively. The 
solid squares mark points for the QSO. The plots show points only for those objects with S/N > 15. 
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Fig. 11 — Plots for the FOR J0238 - 3443 field of the position residuals, TLX and Tiy, as a function of the 
pixel phase, & x and Q y . The panels a, 6, and c are for the epochs 2000, 2001, and 2003, respectively. The 
solid squares mark points for the QSO. The plots show points only for those objects with S/N > 15. 
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RY in the figure indicates the epoch of the measurement. None of the plots shows a trend between RX and 
X or RY and Y. Although this article does not show them, the plots of RX versus Y and RY versus X also 
do not show any trends. Figure [T3l plots weighted mean residuals, (RX) and (RY), in 10 equal- length bins 
in X or Y . The error bar on a point is the usual uncertainty in the weighted mean. The values of (RX) and 
(RY) do not show systematic trends for X and Y greater than approximately 150 pixel. However, the two 
points with X < 150 pixel or Y < 150 pixel show a combination of large error bars due to a small number of 
points contributing to the mean and a deviation from (RX) = pixel or (RY) = pixel due to a mismatch 
between the ePSF and PSF near the edges of the chip. The points with the largest X or Y also tend to show 
a greater than expected deviation, due to the mismatch. The position residuals do not show any systematic 
trends with respect to S/N, as depicted by Figure [Mj The plotted values of the S/N are from the first-epoch 
data. As expected, the scatter in the values of RX and RY increases with decreasing S/N. 

In the standard coordinate system, in which the stars of Fornax are at rest, the QSO moves. Figure [T5l 
depicts such motion for image A. The implied motion from the best- fitting straight lines is (ji x ,fi y ) = 
(-0.0084 ± 0.0019, -0.0102 ± 0.0014) pixel yr -1 . The contribution to the total X 2 from the QSO is 1.47 
for two degrees of freedom. Figure [TBI is analogous to Figure [TBI for image B. Here, the implied motion is 
(fix, fly) = (-0.0043 ± 0.0020, -0.0132 ± 0.0024) pixel yr" 1 and the contribution to the total \ 2 is 1-37 for 
two degrees of freedom. For ease of comparison, the scales of Figures [15] and [16] are identical. 

5.2. Motion of the QSO in the FOR J0240 3438 field 

The number of objects with a measured centroid is 442, 407, and 431 in epochs 2000, 2001, and 2002, 
respectively. Among these, 310 objects are common to the three epochs. The greatest difference in the 
pointing of HST is between the 2000 and 2002 epochs, where the difference is about 2.5 pixel in the X 
direction and about 5.3 pixel in the Y. Because the images used in the analysis were corrected for the effects 
caused by the decreasing charge transfer efficiency, the transformation to the standard coordinate system 
did not include a term that would account for such effects. The choice for the individual % 2 that triggers 
fitting for uniform linear motion is 15. The multiplicative constant that produces a \ 2 of one per degree of 
freedom is 1.100. 

Figure [17] plots RX and RY as a function of position in the standard coordinate system. Although the 
plots do not show any egregious trends, some subtle structures can be discerned. The points in the plot of 
RX versus X for the 2001 epoch (middle row) have a greater scatter than those in the corresponding plots 
for the two other epochs and show a "flare up" around X w 300 pixel. The points in the plot RY versus Y 
also have a greater scatter, which, in contrast, varies little with Y. These observations are more apparent in 
Figure [TS] Only the plot of (RX) versus X for the middle epoch shows a statistically significant deviation 
of (RX) from 0. Here, (RX) for X between 300-400 pixel deviates more than 4 a from 0. Figure [19] shows 
no systematic trends between the position residuals and S/N. 

Figure [201 depicts the motion of the QSO in the standard coordinate system. The implied motion from 
the best-fitting straight lines is (fi x ,fi v ) = (0.0124 ± 0.0030, -0.0002 ± 0.0032) pixel yr" 1 . The contribution 
to the total \ 2 from the QSO is 0.98 for two degrees of freedom. 
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Fig. 12. — Position residuals denned by Equations (10) and (11) in iPiatek et alj (|2005l ) for the objects in 
the FOR J0240 — 3434 field. The epochs increase chronologically from top to bottom. Left: RX vs. X. 
Right: RY vs. Y. The squares correspond to image A of the QSO and the triangles to image B. The arrows 
indicate points beyond the boundaries of the plot. 
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Fig. 13. — Mean position residuals as a function of position for the FOR J0240 — 3434 field. Left: (RX) 
vs. X. Right: (RY) vs. Y. The points are the weighted mean residuals in 10 equal-length bins in X or Y. 
The points are plotted at the mean of the coordinate values in the bin. The numbers in parenthesis indicate 
how many values contributed to the mean. The larger squares and triangles mark those points that include 
contributions to the mean from images A and B of the QSO, respectively. 
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Fig. 14. — Position residuals as a function of S/N for the objects in the FOR J0240 — 3434 field. The epochs 
increase chronologically from top to bottom. Left: RX vs. S/N. Right: RY vs. S/N. The squares correspond 
to image A of the QSO and the triangles to image B. The arrows indicate points beyond the boundaries of 
the plot. 
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Fig. 15. — Location of image A of the QSO as a function of time for the FOR J0240 — 3434 field in the 
standard coordinate system. 
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Fig. 16. — Location of the image B of the QSO as a function of time for the FOR J0240 — 3434 field in the 
standard coordinate system. 
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Fig. 17. — Position residuals for the objects in the FOR J0240— 3438 field. The epochs increase chronolog- 
ically from top to bottom. Left: RX vs. X. Right: RY vs. Y. The squares correspond to the QSO. The 
arrows indicate points beyond the boundaries of the plot. 
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Fig. 18. — Mean position residuals as a function of position for the FOR J0240 — 3438 field. Left: (RX) vs. 
X. Right: (RY) vs. Y. The points arc the weighted mean residuals in 10 equal-length bins in X or Y. The 
points are plotted at the mean of the coordinate values in the bin. The numbers in parenthesis indicate how 
many values contributed to the mean. The larger squares mark those points that include the contribution 
to the mean from the QSO. 
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Fig. 19. — Position residuals as a function of S/N for the objects in the FOR J0240 - 3438 field. The epochs 
increase chronologically from top to bottom. Left: RX vs. S/N. Right: RY vs. S/N. The squares correspond 
to the QSO. The arrows indicate points beyond the boundaries of the plot. 
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Fig. 20. — Location of the QSO as a function of time for the FOR J0240 - 3438 field in the standard 
coordinate system. The vertical axis in each plot has the same scale. 

5.3. Motion of the QSO in the FOR J0238 - 3443 field 

The number of objects with a measured centroid is 365, 380, and 347 in epochs 2000, 2001, and 2003, 
respectively, and 256 are common to all three epochs. The greatest difference in the pointing of HST is 
between the 2001 and 2003 epochs, where the difference is about 2.6 pixel in the X direction and about 
15.6 pixel in the Y . As for the previous field, the images were corrected for the effects caused by the 
decreasing charge transfer efficiency. The choice for the individual x 2 that triggers fitting for uniform linear 
motion is 15. The multiplicative constant that produces a x 2 of one per degree of freedom is 1.125. 

Figures [21] and [22] for the FOR J0238 - 3443 field are analogous to Figures [17] and [18] The first of these 
two figures does not show any alarming trends between the residuals and position. The scatter of the points 
shown in the bottom row of plots — especially in (RX) versus X — is larger than that in the other plots. 
The larger error bars and somewhat larger scatter around zero seen in the plot of average residuals confirm 
this observation. Figure [221 analogous to \T§\ also does not show any systematic trends between the position 
residuals and S/N. 

Figure [24] shows the motion of the QSO in the standard coordinate system. The implied motion from 
the best-fitting straight lines is (ji x ,fi y ) = (0.0103 ± 0.0011, -0.0056 ± 0.0014) pixel yr" 1 . The contribution 
to the total \ 2 from the QSO is 0.33 for two degrees of freedom. 



5.4. Measured Proper Motion 

Table 3 gives the measured proper motion for each field in the equatorial coordinate system and the 
weighted mean. Table 4 gives the proper motion for those objects in the FOR J0240 — 3434 field for which 
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Fig. 21. — Position residuals for the objects in the FOR J0238 — 3443 field. The epochs increase chronolog- 
ically from top to bottom. Left: RX vs. X. Right: RY vs. Y. The squares correspond to the QSO. The 
arrows indicate points beyond the boundaries of the plot. 
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Fig. 22. — Mean position residuals as a function of position for the FOR J0238 — 3443 field. Left: (RX) vs. 
X. Right: (RY) vs. Y. The points arc the weighted mean residuals in 10 equal-length bins in X or Y. The 
points are plotted at the mean of the coordinate values in the bin. The numbers in parenthesis indicate how 
many values contributed to the mean. The larger squares mark those points that include the contribution 
to the mean from the QSO. 
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Fig. 23. — Position residuals as a function of S/N for the objects in the FOR J0238 - 3443 field. The epochs 
increase chronologically from top to bottom. Left: RX vs. S/N. Right: RY vs. S/N. The squares correspond 
to the QSO. The arrows indicate points beyond the boundaries of the plot. 
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Fig. 24. — Location of the QSO as a function of time for the FOR J0238 
coordinate system. The vertical axis in each plot has the same scale. 



3443 field in the standard 



it was measured. Tables 5 and 6 do the same for the FOR J0240 - 3438 and FOR J0238 - 3443 fields, 
respectively. The first two lines in Table 4 and the first lines in Tables 5 and 6 correspond to the QSOs. 
The subsequent lines list objects in order of decreasing S/N. After column (1) listing the ID of an object, 
columns (2) and (3) give the X- and y-coordinates of the object in the earliest image of the first epoch. 
These images are u50j0201r, o5bl04010, and o5bl03010 for Tables 4, 5, and 6, respectively. Column (4) gives 
the average S/N at the first epoch. Columns (5) and (6) give the equatorial components of the measured 
proper motion. Each value is the sum of the measured motion in the standard coordinate system and the 
weighted mean proper motion of Fornax given in the bottom line of Table 3. In other words, the listed values 
are with respect to the QSOs and, thus, the values for the QSOs are listed as 0. The quoted uncertainty is 
the sum in quadrature of the uncertainty in the motion of an object in the standard coordinate system and 
the uncertainty of the weighted mean proper motion of Fornax. Column (7) gives the contribution of the 
object to the total \ 2 - 



The measured proper motion of Forn ax from this article, ( n, 



(47.6±4.6,- 36.0±4.1) mas century' 



replaces the preliminary measurement in iPiatek et al.l (j2002f ). iDinescu et al.l (|2004l ) reported another inde- 
pendent measurement, (/i Q , fig) = (59 ± 16, —15 ± 16) mas century -1 . For these two independent measure- 
ments, the a components differ by 0.68 times the uncertainty of their difference, whereas the d components 
differ by 1.27 times. 



Figure [55] compares the measurement from lDinescu et al.l (|2004) and the four independent measurements 
from this article. A rectangle represents a single measurement: its center, marked with a circle, is the best 
estimate of the pro per motion and the s ides are offset by the 1-er uncertainties. Rectangle 1 represents the 



measurement from IDinescu et al.l (|2004j ). while the other four represent the measurements in this article 
(see the caption). The a components of all the measurements differ from their mean by less than their 
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Fig. 25. — Comparison of five independent measurements of the proper motion of Fornax. The center of each 
rectangle, marked with a circle, is the best estimate of the proper moti on, and the sides are offset by the 1 
a uncertainties. Rectangles 1-5 correspond to the measurements from lDinescu et al.l (|2004r ). and from this 
study for the fields FOR J0240 - 3434A, FOR J0240 - 3434B, FOR J0240 - 3438, and FOR J0238 - 3443, 
respectively. 



uncertainties, whereas the 5 components of some measurements differ by more. For the measurements 
from this study, the largest difference in the S components, for 2 and 3, is 1.55 times its uncertainty. The 
X 2 of our measurement s around their weigh ted mean is 1.3 for \i a and 2.9 for fig, each for 3 degrees of 



freedom. Including the iDinescu et al.l (|2004T ) measurement in the weighted mean alters the values by only 
1-2 mas century" 1 . The \ 2 of the five measurements around this mean is 1.9 for fi a and 4.8 for (x$, each 
for 4 degrees of freedom. We conclude that the agreement among the five independent measurements is 
reasonable, implying that the proper motion of Fornax is reliably known. 



The result in lDinescu et al.l (|2004h is a weighted mean of two measurements, one using g alaxies and the 
other QSOs as standards of rest. The lensed QSO J0240 — 3434 is among the standards, and lDinescu et al 



(|2004l ) give individual values of the proper motion with respect to image A and B of the QSO. They are: 
(a*q ,^s) = (84 ± 52, — 4 ± 64) mas century -1 for image A and (fx a ,fJ,s) = (73 ± 34, — 93 ± 80) mas century" 1 
for image B. Comparing these to our measurements in Table 3 shows an agreement within the uncertainties; 
however, because of the large uncertainties in the values above they do not provide a significant check on 
our results. 
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5-4.1. Effect of Changing Charge Transfer Efficiency 

We have performed a parallel analysis that made no correction for the effects of charge transfer efficiency 
for each of the three fields. The resulting weighted mean fi a is nearly the same as the value in Table 3, but 
the mean /j,s is more negative by 17 mas century" 1 . The \ 2 around these means are 3.0 and 2.9, respectively, 
both for 3 degrees of freedom. These values are larger than those for the results in Table 3, demonstrating 
that correcting for the effects of the changing charge transfer efficiency significantly improves the agreement 
between the individual measurements of the proper motion. 



5.5. Galactic Rest-Frame Proper Motion 

The measured proper motion of the dSph contains contributions from the motion of the LSR and the 
peculiar motion of the Sun. The magnitude of each contribution depends on the location of the dSph in the 
sky. Removing them produces the Galactic rest-frame proper motion — the proper motion measured by a 
hypothetical observer at the location of the Sun but at rest with respect to the Galactic center. Columns 
(2) and (3) of Table 7 give the equatorial components, (^° rf , ^° rf ), of the Galactic-rest-frame proper motion. 
Their derivation assumes: 220 km s _1 for the circular velocity of the LSR; 8.5 kpc for the distance of 
the Sun from the Galac tic center; and (u g ,v q ,Wq) = (-10.00 ± 0.36, 5.25 ± 0.62, 7.17 ± 0.38) km s" 1 



([Dehnen fc Binnevlll998f ) for the peculiar velocity of the Sun, where the components are positive if Uq points 
radially away from the Galactic center, vq is in the direction of rotation of the Galactic disk, and wq points 
in the direction of the north Galactic pole. Columns (4) and (5) give the Galactic rest-frame proper motion 
in the Galactic coordinate system, (/uf rf , //° rt )- Columns (6) - (8) give the II, 9, and Z components of the 
space velocity in a cylindrical coordinate system centered on the dSph. The components are positive if II 
points radially away from the Galactic axis of rotation, points in the direction of rotation of the Galactic 
disk, and Z points in the dir ection of the north Galactic pole. The derivation of the se components assu mes 



138 kpc ( Saviane et al. 2000l ) for the heliocentric distance to and 53.3 ± 0.8 km s 1 ( Walker et al. 2006h for 



the heliocentric radial velocity of Fornax. Columns (9) and (10) give the radial and tangential components 
of the space velocity for an observer at rest at the Galactic center. The component V r is positive if it points 
radially away from the Galactic center. The uncertainties in the listed quantities derive from Monte Carlo 
experiments. The bottom line in Table 7 gives the weighted mean of each listed quantity. Note that even 
though II 2 + 6 2 + Z 2 = V 2 + V 2 for a measurement in a given field, this equality may not hold for the mean 
values. 

The negative V r in Table 7 means that Fornax is moving towards the Milky Way. The negative value 
for 9 shows that the orbit is retrograde. 



6. Orbit and Orbital Elements of Fornax 

Knowing the space velocity of a dSph allows a determination of the orbit of the galaxy for a given form 
of the Galactic potential. This work adopts a Galactic potential that has a contribution from a disk of the 
form ( Miyamoto fc Nagai 1975h 



*,i^ = ; GMdisk (1) 

R 2 + (a + ^Z 2 +b 2 ) 2 
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from a spheroid of the form (|Hernquistl ll990) 



GAL 



spher 



spher 



^GC 



and from a halo of the form 



*halo =^aloM#GC + rf2 )- 



(2) 



(3) 



In the above equations, Rqq is the Galactocentric distance, R is the projection of Rqq onto the plane of 
the Galactic disk, and Z is the distance from the plane of the disk. All ot her quantities in the eq uations are 
adjustable parameters and their values are the same as those adopted by Uohnston et al.l (|1995T ): Af^isk = 



1.0 x 10 11 M , M, 
d = 12.0 kpc. 



spher 



3.4 x 10 10 M 



©' v halo 



128 km s" 



6.5 kpc, b — 0.26 kpc, c = 0.7 kpc, and 



Figure [26] shows the projections of the orbit of Fornax onto the X-Y (top left panel), X-Z (bottom left 
panel), and Y-Z (bottom right panel) Cartesian planes. The orbit results from an integration of the motion 
in the Galactic potential given by equations (H}-©. The integration extends for 6 Gyr backward in time 
and begins at the current location of Fornax with the negative of the space velocity components given in 
the bottom line of columns (6)-(8) of Table 7. The squares in Figure [2U mark the current location of the 
dSph, the stars indicate the center of the Galaxy, and the three small circles mark the points on the orbit 
at which Z = 0, or, in other words, where the orbit crosses the plane of the Galactic disk. The large circle 
is for reference; it has a radius of 30 kpc. In the right-handed coordinate system of Figure [26j the current 
location of the Sun is on the positive X-axis. The figure shows that Fornax is moving toward the Milky 
Way, is closer to apogalacticon than perigalacticon, and that it has a moderately polar orbit with a small 
eccentricity. 

Table 8 tabulates the elements of the orbit of Fornax. The value of the quantity is in column (4) and 
its 95% confidence interval is in column (5). The latter comes from 1000 Monte Carlo experiments, where 
an experiment integrates the orbit using an initial velocity that is generated by randomly choosing the line- 
of-sight velocity and the two components of the measured proper motion from Gaussian distributions whose 
mean and standard deviation are the best estimate of the quantity and its quoted uncertainty, respectively. 
The eccentricity of the orbit is defined as 

e _ (Rg - Rp) (4) 

(R a + Rp) 

The most likely orbit has about a 1.3:1 ratio of apogalacticon to perigalacticon, and the 95% confidence 
interval for the eccentricity allows ratios approximately between 1.2:1 and 2.2:1. The orbit is retrograde and 
only 11° away from bein g polar. The orbital period, 3.2 Gyr, is about 45% longer than that for Sculptor (2.2 



Gyr; Piatek et al.l (|2006f and mo re than a fact o r of tw o longer than those for Carina (1.4 Gvr: iPiatek et al 
42003J)) and Ursa Minor (1.5 Gvr: |Piatek et all 120051 )1. 



7. Discussion 
7.1. Is Fornax a Member of a Stream? 



Lvnden-Belll (| 19821 ) proposed that most, if not all, dSphs have a "tidal origin." In this picture, the 



tidal field of the Galaxy distorts a progenitor galaxy and creates streams of stars and gas which become 
gravitationally unbound to the progenitor galaxy, but move on a similar orbit. A new dSph forms from this 
material — tidal debris — when dissipation of energy makes distinct clumps of the debris gravitationally 
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Fig. 26. — Projections of the orbit of Fornax on the X-Y plane (top left), the X-Z plane (bottom left), and 
the Y-Z plane (bottom right). The origin of the right-handed coordinate system is at the Galactic center, 
which is marked with a star. The Galactic disk is in the X-Y plane and the present location of the Sun is 
on the positive X axis. The squares mark the current location of Fornax at (X, Y, Z) = (39, 48, —126) kpc. 
For reference, the large circle in the X-Y plane has radius of 30 kpc. The three small circles in the X-Y 
projection mark the points at which Fornax passes through the plane of the Galactic disk. The integration 
starts from the present and extends backward in time for 6 Gyr. 
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bound. Consequently, the galaxies and globular clusters that formed from the deb ris move on simila r 
orbits or, in other words, constitute a "stream." In the case of the "FLS" stream ( Lvnden-Bell 1982), 
Fornax is the survi ving progenitor galaxy, whereas Leo I, Leo II, and Sculptor are the "condensate" dSphs. 
Lvnden-Belll rtl976f) prop osed that another great stream of galaxies is associated with the Magellanic Stream 
(jMathewson et al.lll974l ) and contains the Small and Large Magellanic Clouds, Draco, and Ursa Minor. 
The numbe r of memb e r gal axies in the streams increased as new satellite galaxies of the Milky Way were 
discovered. iMaiewskil <| 19941 ) modified the FLS stream by including in it Sextans and renaming it as the 
"FL 2 S 2 " stream. The orbital plane of the FL 2 S 2 stream is only slightly different from the original plane of 
the FLS stream. The idea of streams is verifiable. Given distances to and radial velocities of the members of a 
stream, the expected proper motion s of the members can be calculated. Thus, in a more quantitative analysis, 
Lvnden-Bell fc Lvnden-Belll (| 19951 ) calculated likelihoods for the existence of several streams containing 
dSphs and globular clusters and predicted the proper motions for the member galaxies and globular clusters. 
Fornax could be a member of one of four possible streams: streams la or lb (together with Pal 14, Pal 15, 
and Eridanus), 4a (together with Sextans, Sculptor, and Pal 3), and 4b (together with Sextans and Sculptor). 
Table 9 contains the predicted motions for Fornax; columns (2) and (3) give the a and 5 components of the 
proper motion, whereas columns (4) and (5) give this motion in polar coordinates. For easy comparison, the 
bottom line lists our measurement from Table 3. 



The proper motion of Fornax predicted for stream 4a is in the closest agreement with our measurement, 
though even here the difference is more than 3.2 times the measurement uncertainty in one component. Given 
that there must be some uncertainty in t he predicted motion s, this difference may not be large enough to rule 
out membership in stream 4a. However. iPiatek et al.l (|2006l ) ruled out the membership of Sculptor in stream 
4a, so the reality of the stream hinges on the agreement between the predicted and actual proper motions 
of the remaining three members. Currently, neither Sextans nor Pal 3 has a measured proper motion. 



7.2. Is the Orbit of Fornax in the Kroupa-Theis-Boily Plane? 



Kroupa et al. ( 20051 ) noted that viewing 11 of the innermost dwarf galaxies, satellites of the Milky Way, 



from a point at infinity with I = 167?91 shows that the galaxies lie on nearly the same plane, whose poles 
are at (l,b) = (168°, -16°) and (348°, +16°). The probability is less than 0.5% that such an alignment of 
galaxies occurs by chance if the galaxies have an isotropic distribution of orbits. If the planar alignment is 
to persist in time, the dwarf galaxies must have orbits contained in this plane. Therefore, its persistence can 
be tested if proper motions of the galaxies are known. 

Let the pole of an orbit be in the direction of the orbital angular momentum vector. Thus, the Galactic 
coordinates of the pole in terms of the orbital elements are: 

(£,b) = (ft + 90°, $-90°). (5) 

For prograde orbits, b < 0, and for retrograde, b > 0. Using the best estimates for ft and $ from Table 8 
and uncertainties that are 1-cr errors from Monte Carlo simulations, the orbital pole of Fornax is at (£, b) = 
(163° ± 8°, 13° ± 4°). The poles of the orbit and plane are separated by 29.4° ± 4.1°. The separation is 
larger than the measurement uncertainties, however, the agreemen t betw een the two poles does not need 
to be exact because the plane has a finite thickness. Kroupa et al. ( 2005h find an rms distance from their 



plane of 26.4 kpc for the sample of 11 galaxies — which extends to about twice the distance of Fornax. This 
thickness is consistent with their claim of a ratio of height to radius of 0.15 for the disk. The tilt between 
the orbit of Fornax and the plane predicts that the time-averaged rms distance of Fornax from the plane is 
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47 kpc. This calculation assumes a circular orbit with a radius of 135 kpc, which is the average of the peri- 
anal apogalacticons. We conclude that the orbit of Fornax is inconsistent with the plane. Whether the orbits 
of the dSphs are consistent with the predictions of A-CDM models for the Local Group will be addressed in 
a later article when measurements become available for Draco and Sextans. 



7.3. What Happened to the Gas in Fornax? 



Stetson et all (|l998h and lSaviane et all |200C| ) detected a population of blue and young main-scqucncc 



stars in Fornax. Th e stars are about 10 s years old and located in the central region of the galaxy. Quoting 
Stetson et al.l (|1998l ). the population "has a flattened distribution on the sky and its major axis is offset by 
roughly 30° from the symmetry axis of the galaxy as a whole." T his dete c tion i mplies that Fornax had to 
have HI 10 s years ago. However, a search for HI with the VLA by Young ( 19991 ) produced no detection at 



the limit of 4. 6 x 10 18 cm 2 a t the galaxy center and 7.9 x 10 18 cm 2 at the distance of one core radius from 



the center. iBouchard et al.l (|2006j ) searched for HI near t he dwar f galax ies in the Local Group using the 
Parkes radio telescope. The survey confirmed the result of Youns ( 19991 ) that there is no HI at the center 
of Fornax, but did identify gas offset about 30' north and west-northwest fro m the center. The str ongest 
emission is between position angles of about 325° and 360° (see Figure 2 of IBouchard et alj|2006l) . The 
velocity of the gas is also offset — about 30 km s" 1 below that of the stars. IBouchard et al.l (|2006r ) calculate 



an approximately 10% probability that this gas is a chance superposition of a high-velocity cloud. This 
probability is not low enough to decisively establish a physical relationship between the gas and the galaxy. 
Bouchard et al. ( 20061 ) also note that subtracting the emission from the Milky Way is difficult for these 
observations. If there is no HI associated with Fornax, then why do young stars reside in the central region? 
If the HI is associated, why is it offset in both position and velocity? 

Trying to answer such questions, Youne ( 19991) offers two possibilities: (1) there is gas in the central 
region but it is either ionized or molecular or (2) HI was displaced from the central region in the past 10 8 
years. There are no published results of searches for HII or molecular gas, and so, (1) cannot be ignored or 
excluded until such results exist. In the case of (2), there must be a physical mechanism f or displacing the gas. 
A su pernova-driven wind blowing the gas from the center is one possible mechanism (|Mac Low fc Ferrara 
19991 ). Ram pres sure due to the mot ion of the galaxy throu gh a gaseous halo of the Milky Way is another 



mechanism, e. 



Gallart et al. ( 2001 ) or Maver et al. ( 20061 ). In this case, the gas would trail the dSph and 
be displaced in the direction opposite to the Galactic rest-frame proper motion vector. This vector has a 
position angle of 120°±8.5°, i.e. pointing in the southeast direction, which would predict that the gas should 
be at a position angle of about 300°. However, the strongest HI emission has a position angl e of ab out 340° 



and the elongation of the cloud is not aligned with the proper motion vector. iDinescu et al.l (|2004j) propose 
a variation of the ram-pressure picture in which Fornax cross ed the Magellanic plane — t he plane containing 
the orbit of the LMC as defined by the space velocity from Ivan der Marel et al.l (J2002J) — about 190 Myr 



ago and passed through gas that was removed from the LMC but which still trails this gas-rich galaxy. This 
picture can be tested using our new proper motion for Fornax and a recently-measured proper motion for 
the LMC. 



Kallivavalil et alJ (|2006al ) measured the proper motion of the LMC using an observational technique 
similar to the one employed in this contribution. The measurement is based on HST observations of 21 
QSOs in the direction of the LMC with the High Resolution Camera of the Advanced Cam era for Surveys. 
The QSOs were discovered by their variability in the MACHO survey (jGeha et al.l 120031 ) . The measured 
proper motion of the LMC is (fi a , fis) — (203 ± 8, 44 ± 5) mas century" 1 . 
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Integrating the two orbits backwards in time in the potential given by Eqs. [TH3] shows that Fornax 
crosses the orbital plane of the LMC about 1.4 Gyr ago and it does so at a distance of only 1.9 kpc from 
the orbit of the LMC. However, when Fornax is closest to the LMC orbit, the LMC will not arrive at that 
point for about another 0.8 Gyr. Numerical experiments that integrate the orbits using starting velocities 
with uncertainties determined by the observations indicate that Fornax crosses the plane of the LMC orbit 
between 0.80 Gyr and 2.0 Gyr ago with 95% confidence. Fornax approaches closer than 20 kpc to the orbit 
of the LMC in about 50% of the simulations. The Monte Carlo simulations show that the two galaxies do 
not come closer than 81 kpc at 97.5% confidence. 

For the orbits adopted here, the recent loss of HI from Fornax 190 Myr ago cannot be attributed to the 
ram pressure from gas in the LMC orbital plane. The last crossing of the orbital plane by Fornax occurred 
about 1.4 Gyr ago, which is 1 Gyr before the cessation of star formation in the dSph. The plane crossing 
occurre d 0.6 Gyr befo r e the cessation of star formation with 97.5% confidence. Our conclusion differs from 



Dinescu et all (|2004) because of differences in the adopted orbits for Fornax and, particularly, the 



that of 
LMC. 

Though an interaction between Fornax and gas from the LMC cannot explain the recent cessation of 
star formation in the dSph, the orbits of the two galaxies are likely to approach closely. It is still an open 
question whether this intersection of orbits has affected the star formation history or gas content of Fornax. 
The strongest argument against such an effect is that the closest approach of the two orbits occurs before 
the arrival of the LMC at that point. However, the am ount of gas in the leading and trailing parts of the 
orbit are uncertain. The proper motion for the LMC of iKallivavalil et all (|2006al ) implies a Galactocentric 
space velocity of 380 km s^ 1 , which yields an orbital period — calculated using the method in Section [6] 
- of 4.3 Gyr. This period is almost t hree times longer than that assumed in the canonical models fo r 
the formation of the Magellanic Stream ( Gardiner 1999t Yoshizawa fc Noguchi 2003 ; Conners et al. 2004 ). 
This conflict between the measured motion of the LMC and the models for the formation of the Magellanic 
Stream means that it is premature to rule out the removal of gas from Fornax by interaction wi th the Stream. 



An in teraction with HI leading the LMC could have removed at least some gas from Fornax. IPutman et al 



(|l998l ) identified what appears to be a leading tidal tail in HI maps. However, the likelihood of an interaction 
is decreased because the leading tail contains only about 25% as much gas as t he trailing tail and the gas is 
organized in three disjoint complexes rather than in a continuous distribution (jBriins et al.1 120051 ) . 



7.4. The Effect of the Galactic Tidal Force on the Structure of Fornax 



If Fornax formed as a spherical system and if the Galactic tidal force has had an important effect on its 
struc ture, then the project e d ellipticity would be the result of the stretching of the dSph in its orbital plane 
(e.g., Johnston et al. 1995 : Oh et al. 1995 : Piatek fc Prvor 199ot ) . Thus, a nearly Galactocentric observer 
would see the Galactic rest-frame proper motion vector aligned with the projected major axis. From Table 1, 
the position angle of the major axis is 41° ± 1°, whereas this angle for the weighted mean Galactic rest-frame 
proper motion is 120° ±8.5°. Within the uncertainties, these two angles are almost orthogonal to each other; 
in other words, the true major axis is not in the orbital plane. 

Notwithstanding the disagreement in the two position angles above, the structure of Fornax exhibits 
several f eatures which c o uld h ave arisen as a consequence of tidal interaction with the Milky Way. Most 
recently, Coleman et al. ( 2005h show that: (1) the ellipticity of the isodensity contours increases with in- 
creasing semi-major axis, reaching a maximum at around 35', and then decreases beyond this radius; (2) the 
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position angle of ellipses fitted to the stellar surface density decreases with increasing semi-major axis; (3) 
isodensity contours on the eastern side of the galaxy are closer together than those on the western; and (4) 
the map of the stellar surface density shows stars of the dSph beyond the fitted tidal radius along the minor 
axis. 

The feature which is most suggestive of a tidal interaction is (4) . The direction of the Galactic rest-frame 
proper motion supports this picture; however, other aspects of the orbit do not: the large perigalacticon, 
long orbital period, small eccentr icity, and approaching motion of the dSph do not favor the tidal origin for 



the observed "extratidal" stars. IColeman et al.l (|2005l ) propose that instead of the Galactic tidal force, a 



collision and merger event with another system is responsible for this observed feature. 



7.5. Is There Dark Matter in Fornax? 



Assuming that mass follows light and using the measured velocity dispersion. iMateo et al.l (| 1 9 9 lh derived 
a M/Ly of 12.3 ± 4.5 in solar units for Fornax, noting that M/Lys in the range between 5.3 and 26 are also 
possible given the un certa inties in the f i tted s tructural parameters. The subsequent derivations of M/L by 
Walcher et all |2003h and IWalker et al.l (|2006f) are within this quoted range. Knowing the orbit of a dSph 
permits another approach to constraining its mass. The mass and orbit determine the tidal radius — the 
radius beyond which a star becomes unbound from the dSph due to gravitational interaction with the Milky 
Way. Making the assumption that the tidal radius is equal to the observed limiting radius then pr ovides 
a valu e for the M/L. The value commonly adopted for the limiting radius is the limiting radius of a [King 
(jl966l ) model fitted to the radial surface brightness profile (unfortunately often also called the tidal radius) . 
An exponential — which has no limiting radius — is also a good fit to the surface brightness profiles of 
dSphs, so the identification of the limiting radius determined from a King model with the tidal radius may 
be incorrect. Nevertheless, we calculate the M/L of Fornax by equating the limiting and tidal radii. A 
constraint on the M/L of a dSph with fewer assumptions, also given below, is that it must be large enough 
for the dSph to have survived the tidal interaction with the Milky Way given its measured orbit. 



Equation ((6]) (King 19 62tlOh et al.lll992n gives the tidal radius, rt, as a function of the orbital elements 
of a dSph moving in a logarithmic potential for the Milky Way. 

( (!- e ) 2 M\ 1/3 



V [(1 + e) 2 /2e] ln[(l + e)/(l - e)] + 1 Mq 

In the equation, e is eccentricity of the orbit, a is the semi-major axis (a = (i? a + i? p )/2), M is the mass of the 
dSph, and Mq is the m ass of the Milky Way within a. Equating r t with the observed limiting radius derived 
by fitting a lKin m odel, rk, gives a value for M/ L for a given orbit. Estimates of the limiting radius 

~~ i). If r k = 64', then 10.4% of the 



200 



for Fornax range from 64' ([Coleman et al.ll2005j ) to 98' (jWalcher et al. 
orbits have M/Ly > 5.3 in 1000 Monte Carlo simulations similar to those that determined the uncertainties 
in the orbital elements; 0.2% of the orbits have M/Ly > 26. If r k = 98', then 100% of the orbits have 
M/Ly > 5.3 and 5.3% have M/Ly > 26. These results show that there is only a 5.3% chance that the 
global M/Ly of Fornax is greater than 26 if the larger of the limiting radii is identified with the tidal radius. 
In other words, the orbit of Fornax is consistent with the broad range of M/L values derived from virial 
equilibrium models and the observed limiting radius of the dSph and, thus, these calculations do not provide 
a strong constraint on the amount of dark matter in Fornax. 

The average measured M/Ly for Galactic globular clusters is 2.3 jPrvor fc Mevlanlll993h . Could Fornax 
have a M/Ly equal to this value and survive the effect of the Galactic tidal force until the present on our 
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derived orbit? Numerical simulations of tidal interactions by lOh et al.l (| 19951 ) and iPiatek fc Prvor (|l995l) 
show that the ratio of the limiting radius derived by fitting a theoretical King model l|Kinelll966f ). to the 
tidal radius defined by Equation ^ is a useful indicator of the importance of the Galactic tidal force on the 
structure of a dSph. These simulations show that: if r^Jr t < 1.0, the Galactic tidal force has little effect on 
the structure of the dSph; at r^/rt « 2.0, the effect of the force increases rapidly with increasing rk/r t ; and, 
for rfe/r f « 3.0, the dSph disintegrates in a few orbits. Assuming that M/Ly — 2.3 and r*fc = 64', r^jr t > 2.0 
for 0.3% of the orbits generated in Monte Carlo simulations. If = 98', the fraction is 10.8%. Thus, it 
is very likely that Fornax could have survived for a Hubble time on its current orbit while containing only 
luminous matter. We conclude that our orbit of Fornax and the observed limiting radius do not constrain 
the amount of dark matter in the dSph. 



7.6. A Lower Limit for the Mass of the Milky Way 



The space velocity of Fornax imposes a lower limit on the mass of the Milky Way, which is given by: 

R(V r 2 + V t 2 ) 



M = 



2G 



(7) 



where R is the Galactocentric distance, and V r and Vt are the radial and tangential components of the 
space velocity, respectively, as measured by an observer at the center of the Galaxy. Equation [7] assumes a 
spherically symmetric mass distribution and for the total energy of the satellite galaxy. Setting R = 140 kpc 
and using the values from Table 7 for V r and Vt, M = (6.4 ± 1.9) x 10 n MQ. This estimate does not place 
a strong constraint on the amount of mass within 140 kpc of the center of the Milky Way since the value is 
similar to the 5.4^^ x 10 n M o within R = 50 kpc found bv lSakamoto et al.l (|2003h . 



8. Summary 

This contribution continues a series of articles, each reporting the proper motion of a dS ph; here, we 



presen t a measurement for Fornax which replaces the value based on only part of the data in IPiatek et al 



(|2002l ). The data come from imaging with the HST and STIS or PC2. Using the measured proper motion, 
we integrate the motion of the dSph to obtain its orbit and then discuss the proposed membership of Fornax 
in streams of satellites, the possible alignment of orbits of the Galactic satellites in a single plane, the 
association between HI and the dSph, the proposed interaction between the LMC and Fornax, the tidal 
interaction of the dSph with the Milky Way, the dark matter content of Fornax, and a lower limit on the 
mass of the Milky Way. The major findings are itemized below. 

1. The weighted mean of the four independent measurements of the proper motion is ([i ai fis) = 
(47.6 ± 4.6, —36.0 ± 4.1) mas century -1 in the equatorial coordinate system for an observer at and moving 
with the Sun. 

2. The Galactic rest-frame proper motion is that measured by an observer at the location of the 
Sun and at rest with respect to the Galactic center. In this system, the proper motion is (Ma rf ) M? 1 ') = 
(24.4 ± 4.6, -14.3 ± 4.1) mas century -1 . 

3. The weighted means of the radial and tangential components of the space velocity with respect 
to a stationary observer at the Galactic center are V r = —31.8 ± 1.7 km s -1 and Vt = 196 ± 29 km s _1 , 
respectively. The negative sign of V r means that Fornax is approaching the Milky Way. 
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4. The orbit resulting from the integration of the best estimate of the motion in a realistic potential 
of the Milky Way shows that Fornax is approaching a perigalacticon of 118 kpc on an approximately polar 
orbit with an eccentricity of 0.13. The apogalacticon of the orbit is 152 kpc and the orbital period is 3.2 Gyr. 



5. Fornax is not a member of s tream la, lb, or 4b o f lLvnden-Bell fc Lvnden-Belll (jl995l ). though it could 
be a member of stream 4a. Since iPiatek et al.l (|2005f ) ruled out the membership of the Sculptor dSph in 
stream 4a, the reality of this stream depends on whether the proper motions of the other proposed members, 
which still need to be measured, agree with the predictions. 

6. The orbit of Fornax is not in the Kroupa-Theis-Boily plane, which contains eleven of the Galactic 
satellite galaxies. 



7. Our findings do not support the hypothesis of iDinescu et alJ (|2004J l that Fornax lost its HI about 
200 Myr ago as it pas sed through the s t ream o f gas trailing the LMC in its orbit. Our proper motion 
for Fornax and that of iKallivavalil et al.l |2006aT ) for the LMC show that Fornax approaches within about 
1.9 kpc of the orbit of the LMC, but 1.4 Gyr ago and 0.8 Gyr before the LMC arrives at that point. An 
interaction 1.4 Gyr ago does not explain why star formation in Fornax continued until 200 Myr ago, though 
uncertainties in modeling the Magellanic Stream keep this mechanism for the loss of HI from Fornax from 
being definitively ruled out. 
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Table 1. Fornax at a Glance 



Quantity 


Value 


Reference 


(1) ' 


(2) 


(3) 


Right Ascension, a (J2000.0) 


02 : 39 : 53.1 


van den Bergh (1999 ) 


Declination, 5 (J2000.0) 


-34 : 30 : 16.0 


van den Bergh (1999) 


Galactic longitude, I 


237.245° 




Galactic latitude, b 


-65.663° 




Heliocentric distance 


138 ± 8 kpc 


Saviane et al. (2000) 


Luminosity, Ly 


(1.4 ±0.4) x 10 7 Lq 


Irwin & Hatzidimitriou (1995) 


Ellipticity, e 


0.30 ±0.01 


// 


Position angle 


41° ± 1° 


// 


Core radius 


13.7' ± 1.2' 


// 


Tidal radius 


71.1' ±4.0' 


// 


Heliocentric radial velocity 


53.3 ±0.8 km s" 1 


Walker et al. (2006) 



Table 2. Fornax: Information about Fields and Images 







R.A. 


Deck 


Date 






rp a 

J- exp 


Field 




(J2000.0) 


(J2000.0) 


yyyy — mm — dd 


Detector 


Filter 


(s) 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


FOR J0240 - 


- 3434 


02 40 07.70 


-34 34 20.01 


1999 - 03 - 10 


PC2 


F606W 


18 x 160 










2001 - 03 - 08 






16 x 160 










2003 - 03 - 08 






16 x 160 


FOR J0240 - 


-3438 


02 40 38.70 


-34 38 58.00 


2000- 01 -31 


STIS 


50CCD 


24 x 192 










2001 - 01 - 25 






24 x 192 










2002 - 01 - 29 






24 x 192 


FOR J0238 - 


-3443 


02 38 43.80 


-34 43 53.00 


2000 - 03 - 08 


STIS 


50CCD 


24 x 193 










2001 - 03 - 08 






24 x 192 










2003 - 03 - 09 






24 x 192 


a Number of 


images times the mean 


exposure time. 


The actual exposure 


times of individual 


images may 



differ by a few percent from the mean due to constraints imposed by the orbit of HST. 
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Table 3. Measured Proper Motion of Fornax 







Ma 


M<5 


Field 




(mas century x ) 


(1) 




(2) 


(3) 


FOR J0240- 


3434A 


54.1 ±8.5 


-27.5 ±7.1 


FOR J0240- 


3434B 


42.4 ±9.6 


-47.7 ± 10.9 


FOR J0240- 


3438 


53.6 ± 15.8 


-32.5 ± 16.2 


FOR J0238 - 


3443 


44.6 ±7.2 


-39.1 ±6.0 


Weighted mean 


47.6 ±4.6 


-36.0 ±4.1 



Table 4. Measured Proper Motions For Objects in the FOR J0240 - 3434 Field 





X 


Y 




Ma 


M<5 




ID 


(pixels) 


(pixels) 


S/N 


(mas century -1 ) 


(mas century -1 ) 


x 2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1 


458 


461 


163 


0± 10 


0±8 


1.47 


2 


562 


540 


122 


0± 11 


0± 11 


1.37 


3 


162 


463 


60 


6±9 


-28 ± 12 


1.44 



Table 5. Measured Proper Motions For Objects in the FOR J0240 - 34348 Field 





X 


Y 




Ma 


M<5 




ID 


(pixels) 


(pixels) 


S/N 


(mas century -1 ) 


(mas century -1 ) 


x 2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1 


480 


520 


218 


0± 16 


0± 16 


0.98 


2 


781 


960 


24 


116 ±34 


24 ±27 


0.81 


3 


126 


358 


15 


442 ± 54 


-317 ±47 


0.31 
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Tabic 6. Measured Proper Motions For Objects in the FOR J0238 - 3443 Field 





X 


Y 




Ma 


M<5 




ID 


(pixels) 


(pixels) 


S/N 


(mas century -1 ) 


(mas century -1 ) 


x 2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1 


499 


541 


167 


0±8 


0± 7 


0.33 


2 


739 


874 


135 


48 ± 10 


-80 ± 12 


2.65 


3 


905 


278 


111 


4±8 


-29 ±9 


0.48 


4 


719 


833 


104 


28 ±9 


-63 ± 12 


1.24 


5 


955 


996 


27 


52 ±21 


-143 ±22 


4.24 


6 


708 


995 


10 


113 ±61 


-226 ±40 


3.15 



Tabic 7. Galactic-Rest-Frame Proper Motion and Space Velocity of Fornax 







Grf 

fa 


Grf 

Ms 


Grf 
Mi 


Grf 

Mb 


n 


e 


Z 




v t 


Field 




(mas cent x ) 


(mas cent 1 ) 


(km s- 1 ) 


(km s^ 1 ) 


(km s- 1 ) 


(km s" 1 ) 


(km s- 1 ) 


(1) 




(2) 


(3) 


(4) 


(5) 


v°) 


{') 








FOR J0240- 


3434,4 


30.9 ±8.5 


-5.8 ±7.1 


0.2 ±7.2 


31.4 ±8.5 


172 ±51 


-21 ±47 


116 ±23 


-27.7±3.1 


206 ± 54 


FOR J0240- 


3434B 


19.2 ±9.6 


-26.0 ± 10.8 


22.1 ±10.8 


23.5 ±9.6 


109 ±57 


-158 ±71 


94 ±26 


-36.6 ±4.1 


211 ±66 


FOR J0240 - 


3438 


30.4 ± 15.8 


-10.8 ± 16.2 


5.2 ± 16.2 


31.8 ± 15.8 


170 ± 94 


-54 ± 110 


118 ±43 


-29.8 ±6.2 


212 ±96 


FOR J0238 - 


3443 


21.4 ±7.2 


-17.3 ±6.0 


12.9 ±6.1 


24.3 ±7.2 


121 ±43 


-100 ±39 


96 ± 19 


-33.1 ±2.7 


180 ± 44 


Weighted mean 


24.4 ±4.6 


-14.3 ±4.1 


9.6 ±4.1 


26.8 ±4.6 


137 ±27 


-80 ±27 


103 ± 12 


-31.8 ± 1.7 


196 ±29 



I 

to 

I 
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Tabic 8. Orbital elements of Fornax 



Quantity 


Symbol 


Unit 


Value 


95% Conf. Interv. 


(1) 


'(2) 


(3) 


(4) 


(5) 


Perigalacticon 


R p 


kpc 


118 


(66, 137) 


Apogalacticon 


Ra 


kpc 


152 


(144,242) 


Eccentricity 


e 




0.13 


(0.11,0.38) 


Period 


T 


Gyr 


3.2 


(2.5,4.6) 


Inclination 


<$> 


deg 


101 


(94, 107) 


Longitude 


n 


deg 


73 


(58,90) 



Table 9. Predicted Proper Motion of Fornax 





Ma 


M<5 


ImI 


PA 


Stream No. 


(mas century -1 ) 


(mas century -1 ) 


(mas century" 1 ) 


(deg) 


(1) 


(2) 


(3) 


(4) 


(5) 


la 


30 


-15 


34 


117 


lb 


15 


-27 


31 


151 


4a 


38 


-23 


44 


121 


4b 


6.3 


-19 


20 


162 


Our Result 


47.6 ±4.6 


-36.0 ±4.1 


60.0 ±4.4 


127 ±4.1 



